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ABSTRACT 

The prompt emission of the long, smooth, and single-pulsed gamma-ray burst, GRB 141028A, is 
analyzed under the guise of an external shock model. First, we fit the 7-ray spectrum with a two- 
component photon model, namely synchrotron+blackbody, and then fit the recovered evolution of the 
synchrotron vF v peak to an analytic model derived considering the emission of a relativistic blast- 
wave expanding into an external medium. The prediction of the model for the vF v peak evolution 
matches well with the observations. We observe the blast-wave transitioning into the deceleration 
phase. Further we assume the expansion of the blast-wave to be nearly adiabatic, motivated by 
the low magnetic field deduced from the observations. This allows us to recover within an order of 
magnitude the flux density at the vF v peak, which is remarkable considering the simplicity of the 
analytic model. Across all wavelengths, synchrotron emission from a single forward shock provides 
a sufficient solution for the observations. Under this scenario we argue that the distinction between 
prompt and afterglow emission is superfluous as both early and late time emission emanate from the 
same source. While the external shock model is clearly not a universal solution, this analysis opens 
the possibility that at least some fraction of GRBs can be explained with an external shock origin of 
their prompt phase. 

Subject headings: gamma-ray burst: individual (141028A) - radiation mechanisms: non-thermal 
radiation mechanisms: thermal 


1. INTRODUCTION 

Identifying the origin of the dynamical evolution of 
gamma-ray burst (GRB) outflows is an unsolved issue, 
critical to the understanding of both the energetics and 
spectra of these events. One idea is that the emission 
is the result of synchrotron radiation from an external 
forward shock pr opagating into the external circumburst 


medium (CBM) (]Cavallo & Rees|1978 Rees & Meszaros 

1992; |Meszaros & Kees||1993 

(Jhiang &; Dermer 1999; 

Dermer et al. 1999; Dermer & 

Vlitman 1999). This mech- 

anism should 

produce smooth 7-ray pulses with dura- 


tions on the order of a few seconds for typical GRB pa¬ 
rameters. However, the short-time variability (on the 
order of a few milliseconds) of many GRB light curves 
ruled out this postulate as a universal mechanism (e.g. 


Sari fe Piran||1997l |Kobayashi et al.||l997| | Walker et cd. 
2000| and gave favor to several alternative hypotheses 
including models that consist of rapid internal shocks in 


an unsteady outflow 

Rees & Meszaros 1994) or mag- 

netic reconnection (e.g., Spruit et al. 2001fDrenkhahn 

& Spruit 

2002 

/hang & Yan 2011) to account tor the 


observed emission. 
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Still, there do exist long and temporally smooth GRBs 


onds (see 

Golkhou & Butler 

2014 

Golkhou et al. 

2015) 

that do not violate the varia 

bility constraints ol t 

ne ex- 


ternal shock model and can be tested via their spec¬ 
tral evolution as to whether they conform to the well- 
established predictions made by this model. The sim¬ 
plicity of the model affords it the ability to be tested 
both spectrally and temporally, a feature unique to the 
external shock model. The dynamics and spectra of 


the internal shock mo del have been simulated (Daigne 


& Mochkovitch 1998), but the dynamics rely on as- 


sumed and degenerate configurations from variations in 
the wind (e.g. the radial distribution of the Lorentz fac¬ 
tors), forbidding the formulation of unique predictions 
that can be identified in the data. Therefore, it is cur¬ 
rently impossible to test the internal shock model in the 
manner presented here without severe degeneracies. Ad¬ 
ditionally, the internal shock model has trouble efficiently 
converting the internal kinetic energy of a GRB into ra¬ 
diation, which is challenging w hen trying to explain t he 
extreme luminosities observed (Kobayas hi et al. [1997 1. 

Herein, we analyze the bright, long, single-pulsed 
GRB 141028A and find several clues for an ex¬ 
ternal shock origin of its emission. We fit the 
GRB’s time-resolved spectr a with a slow-cooled syn- 
chrotron+blackbody model (Burgess et al. 2014b) and 
examine the evolution of the spectra. The evolution of 
the synchrotron vF v peak (F n ) is fit with_an analytic 
physical model predicted by Dermer et al. (1999). From 
this fit, we obtain physical parameters such as the coast¬ 
ing Lorentz factor and CBM radial profile which can then 
be used to predict how the flux of the prompt emission 
should evolve. Comparing these predictions to the data 
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enables us to test the validity of the model in several 
ways. 

The article is organized in the following manner. In 
Section [2j we introduce the formalism of the external 
shock model to derive a function for E p (t) to fit to data. 
In Section [3j the observations and spectral analysis are 
introduced. Sectionffldetails our application of the exter¬ 
nal shock model to the data. The parameters resulting 
from the analysis are used to make further predictions 
regarding the flux evolution of the prompt phase of the 
GRB. From there, we analyze the photospheric emission. 
In Sections |5.1| and [572] we discuss the compatibility of 
our results with the high and low-energy late-time obser¬ 
vations. 

2. THE EXTERNAL SHOCK MODEL 

The external shock model is built upon the blast-wave 
evolution derived in Blandford & McKee (1976) which 
tracks the evolution of a relativistically expanding fire¬ 
ball into an external medium. The equations can be ap¬ 
plied to GRBs by assuming some fraction of the electrons 
in the shocked external medium is accelerated to high en¬ 
ergies by the shock wave and radiates a fraction of t he 
kinetic energy away via synchrotron radiation (Cavallo 


& Rees 1978 Rees & Meszaros 1992). In this work, we 
use th e analytic formalism developed in |Dermer et al. 
(1999) to fit the spectral evolution of the emission. We 


briefly re view the main equat i ons required and refer the 
reader to Dermer et al. (19991; Chiang & Dermer (1999) 
for more details on the model. 

The blast-wave is assumed to expand into an external 
circumburst medium (CBM) with a radial density evolu¬ 
tion modeled as a power law, 


n(x) = n d x 11 cm 


-3 


( 1 ) 


where no is the initial density and 77 describes the radial 
morphology of the CBM such that rj = 0 is a constant 
density and 77 = 2 describes a stellar wind. The dimen¬ 
sionless radial coordinate is x = r /r d where, following the 
convention that a quantity w = w n 10 n , 


r d = 5.4 ■ 10 


16 


(1 — t ?/3)£’o i 54 


n 0,2^n 


0,2 


73 


cm 


( 2 ) 


is the radius at which the blast has swept up a significant 
r 1 - 1 ) to begin decelera ting (jRees fe 


Meszaros 

1992) 

McKee (l 

976) 


1992 ). Following the solution of Blandford & 


consisting of a coasting phase followed by a deceleration 
phase: 


^ ^ r°a;-9 


x < 1 
1 < x 


( 3 ) 


where To is the coasting Lorentz factor and g is the radia¬ 
tive regime index. For a constant density (77 = 0) CBM, 
g = 3 , 3 /2 indicate the fully radiative and no n-radiative 
(adiabatic) e xpansion regimes respectively (Blandford & 
McKee|1976 ). A fraction of the blast energy is dissipated 
in the shock and accelerates electrons to high energies, 
which subsequently radiate this away via s ynchrotron ra- 
diation. Following the parameterization of|Dermer et al. 


(1999), the temporal evolution of the vF v peak energy, 
this synchrotron radiation can be modeled as 


F n , o 


where 


E p (t) — £0 


£a = 15 


F(a) 


1 4 


1 1 / 2 n T4 

1 ' 0 , 2 < 1 - 3 L 0,2 

l + Z 


"f* keV 


keV 


( 4 ) 


( 5 ) 


is the observed E p at the observed deceleration time, 

r d 


U = ^~(l + z). 

1 0 C 


( 6 ) 


Here, c is the speed of light and q parameterizes the 
magnetic field and shock acceleration microphysics such 
that q = [cb ( rs / 4 )] e e- Here 6 b and e e are the magnetic 
and electron equipartition factors and r s ~ 4 is the shock 
compression ratio. For values expected in GRBs and to 
coincide with our choice of slow-cooling synchrotron to 
fit the spectra, if we take £b — 10 -4 and e e < 0.5 then 
we have values of q ~ 10 -3 . While these values of es 
are lower than the typically assumed values of ~ 10 ~ 2 , 
recent studies (Lemoine et al.|2013 Santana et al.||2014| 
find that the values can be much tower and we therefore 
follow these works. For the analysis, we will allow values 
6 b < 10 -2 to cover the range of typically assumed values. 

Additionally, we can write the measured vF v peak lu¬ 
minosity as 


P p (t) = n n 


r 2-r) 


0 < X < 1 


c 2-^-4ff 1 < x < r 


Vs 


where 


(2g - 3 + ^)m p c 3 T^n 0 A 0 _ ± 

Ho cx - - P -erg s . 

2 g(i + zy 


( 7 ) 


( 8 ) 


Here, m p is the proton mass and A 0 is the blast-wave area 
factor. The prop ortionality i n Equ ation [ 8 ] comes from 
the fact that in Dermer et al. (1999) an empirical photon 
model is used to create light curves and here we will 
be using synchrotron emission directly to fit the time- 
resolved spectra of GRB 1^1028A (see Section [3|. 

With these equations, the external shock model can 
predict several observable quantities, namely, the evolu¬ 


tion of E p and P p which we test below in Section |4.2 


The effects on the E p evolution of different parameters 
is shown in Figure [TJ The parameters have the following 
dominant effects: 

• r 0 : sets the break time 

• 77 : sets the slope of early E p decay 

• g: sets the slope of late E p decay. 

These parameters, along with q which acts as a normal¬ 
ization, will be the free parameters that will be deter¬ 
mined by fitting the observed evolution of E p to Equa¬ 
tion [4j In reality, t he range o f val ues for g are determined 
by the value of 77 (|Dermer et al.||1999|. Since we cannot 
know a priori the value of 77 , we treat g as an independent 
parameter. It is important to note that this analytic ap¬ 
proximation to the blast-wave evolution fails to capture 
some aspects that the full numerical solution exhibits 
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Figure 1. Demonstrating the effect of the physical parameters on 
the evolution of E p with time. It is easy to see that To has the 
strongest effect on the timescale of the burst, while 77 affects the 
early evolution of E p and g the late evolution. 


(Chiang & Dermer 1999). Most notable is the smooth 
transition trom coasting to decelerating at the deceler¬ 
ation radius (rd). The formalism we adopt will serve 
as a proof of concept that will be further applied to a 
larger sample and improved upon with a full numerical 
treatment. 


3. OBSERVATIONS OF GRB 141028A 
3.1. Data Acquisition 

Rapid variability poses a problem for the external 
shock model; therefore, bright, long, single-pulsed GRBs 
provide the most viable candidates for being produced 
by external shocks. GRB 141028A is an example of this 
class of GRB (see Figure [ 2 ]). GRB 141028A was discov- 
ered by the Fermi Gamma-ray B urst Monitor (GBM) 
( Roberts||2014| [Meegan et al.||2009 ) and the bright GRB 
triggered an autonomous repomt of the Fermi s pacecraft 
to o p timize the Large Area Telescope (LAT) (Atwood 
et al.||2009|) for follow-up observations. In grou nd anal- 
L AT al so localized and detected the GRB (Bissaldi 
et ah] 20141 until 10 3 s after the trigger. The lo¬ 
calization led to a target of opportunity observation by 
the narrow-field Swift instruments (the X-ray Telescope, 
XRT; and the Ultra violet Optical Telescope, UVOT) 
(Gehrels et al. 2004). The X-ray an d optical aft erglow 
was detected by XBl' (|Kennea||2014|) and UVOT (1 Siegel 
& Pagani 2014), and subsequently by many different 
ground-based facilities, including a measurement of a 
redshift of z=2.332 wit h the Very Large Telescope/X- 
Shooter instrument (Malesani 20141. We collected op- 
tical/NIR photometry from GCJN circulars, an d con¬ 
struc ted a multi-band SED using GROND data (Kann 
20141. and light curves using r’ and i’ filters (GROND, 
TUTTI II : iTro. a| pRl]), P60: |Cenko| pH4|), LCOGTN: 
Kopac| (|2014|)j. The XRT light curve and spectra were 
obtained fro m the XRT Team Repository (Evans et al. 

We triggered a pre-approved late-time 


2007 


2009). 


Chandra target of opportunity observation to constrain 
the properties of the break in the X-ray afterglow hinted 
at by XRT. The 40 ks Chandra observation was ana¬ 
lyzed with CIAO v4.6, yielding a faint detection which 
was converted to flux using the fit to the XRT spectrum. 
The X-ray and optical afterglow observations began at 
~ 10 4 s after the GBM trigger. It is therefore impossible 


to observe the continuous evolution of the flux from the 
prompt to afterglow phase at low-energy. 



10 15 20 25 

Time [s] 


Figure 2. The luminosity lightcurve of GRB 141028A consisting 
of the synchrotron ( purple ) and blackbody ( red ) components. 


3.2. Spectral and Temporal Analysis 

For the prompt emission, we use GBM time-tagged 
event (TTE) data and Fermi LAT low-energy (LLE) 
data to perform spectral analysis. The LLE technique 
is an analysis method designed to study bright transient 
phenomena in the 30 MeV - 1 GeV energy ran ge, and 


was successfully applied to Fer mi- LAT GRBs (Acker- 
mann et al.|2013|) and solar flares (Ackermann et aF 2012 
Ajello et al.|2014 ). The idea behind LLE is to maximize 
the effective area below ~ 1 GeV by relaxing the stan¬ 
dard analysis requirement on background rejection. The 
first five seconds after the trigger time are excluded be¬ 
cause the count rate was too low to constrain a spectral 
model. 


Using the method of Gao et al. 


(2012), we calculated 


the variability components of GRB 141028A during the 
prompt emission. We find the dominant component to 
be Gar = 24.7/(1 + z) s and an insignificant fast compo¬ 
nent with Gar = 3.8/(1 + z ) s. These values place the 
burst safely within the range of what can be expected 
by an external shock. We use time bins of 1 s to be 
sure to bin below the fast variability component. To fit 
the time-resolved spectra, we employ a two-component 
model consisting of synchrotron emission from an incom¬ 
pletely coo led electron distribution (s low-cooled) and a 
blackbody ([Burgess et al.|2011[|2014b|) (see Figure [3]). In 
Burgess et al.| (|2U14b|), it was shown that single pulsed 
GRBs have spectra that are compatible with slow-cooled 
synchrotron sometimes with and sometimes without the 
addition of a blackbody. This motivates our choice of 
the synchrotron+blackbody photon model for fitting the 
spectra of the GRB 141028A. We note that not all GRB 
spectra are compatible with this photon model. There 
exist GRBs such as GRB 090902B with a clear domi¬ 
nant photospheric co mponent originating from subpho - 
tospheric dissipation (Pe’er et al. 2005 Ryde et al.|2010 ). 

The synchrotron model implemented is fully physical 
in the sense that the spectral shape comes from the syn¬ 
chrotron emissivity and electron distribution alone. The 
spectral fitting should n ot be confused w ith fitting the 
empirical Band function (Band et al.Jl993) with its vari- 
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(a) 


(b) 


Figure 3. The time-resolved vFv spectra in 1 s intervals of the synchrotron (5-27 s) (a) and blackbody (5-15 s) (b) components from the 
GBM+LLE fits. The synchrotron spectrum evolves in time from purple to white while the blackbody evolves from red to yellow. 
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able low-energy index (a). The blackbody is only in¬ 
cluded in time bins that statistically require it via a likeli¬ 
hood ratio test. We find that this statistical requirement 
produces a continuous presence of the blackbody for the 
first part of the observation. The data are well fit by the 
two-component model (see Tabled] for a summary of the 
spectral fits) and we observe two key features that moti¬ 
vate our investigation of the GRB under the framework of 
the external shock model: E p evolves as a broken power 
law in time and the blackbody is only significant for a 
duration of ~10 seconds com pared to the total prompt 
duration of ~25 s (see Section 4.41. 

We calculate the total k-corrected energy (|Hogg et al. 
2002) in the synchrotron component in the 30 keV-300 
MeV interval by summation over each time bin: E SYNW = 
4ndl -F SYNCH Ati ~ 9.1 • 10 53 erg. The total isotropic en¬ 
ergy of the burst is estimated to be several times larger. 
In the following, we take the total isotropic energy (ki¬ 
netic + radiative) of the blast-wave to be 2 £j so = 10 55 erg 
while noting that this is an extremely high value. Never¬ 
theless, the GRB is extremely bright and no mechanism 
is known that is efficient enough to convert the entire 
rest mass of the progenitor to radiation. 

For the late time GeV emission, we performed an un¬ 
binned likelihood analysis of the LAT data to recover 
the energy flux (Fg) and photon index (y p h) of the 
emission (see Table ^1) with the gtlike program dis¬ 
tributed with the Fermi ScienceTool^El We selected 
P7REP_S0URCE_V15 photon events from a 15° circular re¬ 
gion centered at the Swift XRT position (R.A. =322?60, 
Dec.=— 0?23, J2000) and within 105° from the local 
zenith (to reduce contamination from the Earth limb). 
Events with measured energy from 100 MeV to 10 GeV 
are included in our analysis (the highest energy event as¬ 
sociated with this GRB has an energy of 3.8 GeV and ar¬ 
rives 157.5 seconds after the GBM trigger time). Further 
details on the LAT analysis are discussed in Appendix [A] 

4. EXTERNAL SHOCK ANALYSIS 
4.1. E p Evolution 

After performing spectral fits to the data we can test 
the external shock model by fitting the evolution of the 
recovered synchrotron E p with Equation [4| The fit is 
performed with a B ayesi an analysis tool built upon the 
MULTINEST (|Feroz et al. 20091 software which allows us 
to fully explore the correlated parameter space of the 
model. The free parameters in the fit are To, r], g , and 
g_ 3 to which we assign flat priors that are consistent 
with physical expectations (T 0 £ {10,1500}, g £ {0,2} 
g £ {0,3} , q -3 £ {0,10}). We also note that it is not 
possible to fit for the density no, which we set at no = 
{ 1 , 10 , 100 } cm ~ 3 with different values mainly affecting 
the value of To recovered. Calculations in the text and 
figures assume no = 100 cm -3 but are complemented 
with calculations at the other values in associated tables. 
Figure [4] shows the fit of the E p evolution and the best-fit 
parameters are detailed in Table [3j 

We find that g ~ 1.3, which is less than what is ex¬ 
pected if the blast-wave decelerates adiabatically (g = 
1.5). This suggests that the blast-wave is still transition¬ 
ing to the asymptotic limits of Equation [3j To test this 

8 We used version 09-34-02 available from the Fermi Science 
Support Center http://fermi.gsfc.nasa.gov/ssc/ 



Time [s] 


Figure 4. Bayesian fit of the E p evolution of GRB I) 1028A. The 
maximum-likelihood point is indicated by the blue curve and pos¬ 
terior samples are in orange. The fitted parameters are indicated 
in Table [3j 

assumption, we numerically solved for the evolution of T 
with radius via the equations of energy and momentum 
conservation and examined the transition phase in Fig¬ 
ure [5j There is clearly a region that corresponds to our 
recovered value of g ~ 1.3. Late (hundreds of seconds 
after the trigger) time observations smoothly connect¬ 
ing the prompt and afterglow emission would allow us 
to measure the asymptotic value of g. Yet, with the re¬ 
covered fit parameters, we can calculate several physical 
properties of the outflow including r<j = 3.4-10 16 cm and 
the resulting predicted P p (t) evolution though we will 
assume g ~ 1.6 to calculate fluxes as this is the value 


reached asymp t otical ly in the near adiabatic case Chi- 


ang & Dermer (19991. These parameter values are used 


in the following sections to gain more insight about the 
GRB and its surrounding environment. 



Figure 5. The blast-wave Lorentz factor evolution from our nu¬ 
merical simulation of the external shock model (blue dashed line). 
We illustrate that we are most likely finding the blast-wave still 
evolving in the transition phase, which accounts for the shallow 
value of g found in the fit of E p . 


4.2. Peak Flux - E p plane 

Using the recovered parameters from the fit, we can 
use Equations [4] and |7| to predict the evolution in the P p - 
Ep plane. InFigure [6j the data from GRB 141028A 
is plotted with the predicted curve from the external 
shock model deri ved from the E p fits. Noting the dis¬ 
cussion in Section 4.1 negative fluxes would be obtained 
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Table 1 

Results of the time-resolved spectral analysis of GBM+LLE data (10 keV - 300 MeV). 


Tstart [s] 

Tatop [s] 

F|7 N0H [erg s 

“ 1 ] b 

E p [keV] c 

pd 

F™ [erg s- 

- 1 ] 

kT [keV] 

6.0 

7.0 

(1.79 ±0.34)- 

10 ~® 

1109 ±336 

10 .00 a 

(1.61 ± 3.43) ■ 

10" 8 

24.55 ± 26.48 

7.0 

8.0 

(3.19 ±0.54)- 

10 “® 

1378 ± 472 

5.74 ± 1.21 

(1.04 ± 0.54) ■ 

io - 7 

25.12 ±4.50 

8.0 

9.0 

(4.99 ±0.59)- 

10 “® 

1383 ±310 

5.21 ± 0.58 

(2.76 ±0.59) ■ 

10- 7 

40.69 ± 5.40 

9.0 

10.0 

(7.25 ± 0.68) ■ 

10 “® 

1504 ± 310 

4.59 ±0.30 

(1.91 ± 0.68) • 

10- 7 

29.98 ±3.94 

10.0 

11.0 

(6.38 ±0.58)- 

10 “® 

738 ±108 

4.11 ± 0.18 

(4.38 ±0.58) • 

10- 7 

60.01 ± 9.61 

11.0 

12.0 

(5.40 ± 0.42) ■ 

10 ~® 

572 ± 79 

4.18 ±0.21 

(3.21 ± 0.42) • 

10- 7 

35.78 ± 4.54 

12.0 

13.0 

(6.26 ± 0.60) ■ 

10" 6 

498 ± 71 

4.09 ±0.20 

(6.96 ±0.60) ■ 

10- 7 

73.11 ± 10.11 

13.0 

14.0 

(6.56 ± 0.42) ■ 

10 -« 

520 ± 55 

4.23 ±0.18 

(3.12 ±0.42) ■ 

10- 7 

37.28 ± 5.62 

14.0 

15.0 

(4.48 ± 0.31) ■ 

10- 6 

310 ± 46 

3.67 ± 0.13 

(1.33 ±0.31) ■ 

10- 7 

21.56 ± 3.78 

15.0 

16.0 

(5.14 ±0.35)- 

10- 6 

288 ± 34 

3.65 ± 0.12 

(2.13 ±0.35) ■ 

10- 7 

31.06 ± 5.68 

16.0 

17.0 

(2.90 ±0.23)- 

10-6 

168 ± 16 

3.69 ±0.16 




17.0 

18.0 

(2.48 ±0.21)- 

10 "® 

175 ± 21 

3.52 ± 0.14 




18.0 

19.0 

(1.75 ±0.20)- 

10 “® 

147 ± 20 

3.68 ±0.22 




19.0 

20.0 

(2.03 ±0.19)- 

10 “® 

123 ± 18 

3.24 ±0.11 




20.0 

21.0 

(1.77 ± 0.18) ■ 

10 “® 

122 ± 18 

3.33 ±0.13 




21.0 

22.0 

(1.42 ±0.19)- 

10 “® 

148 ± 23 

3.77 ± 0.28 




22.0 

23.0 

(1.30 ±0.17)- 

10" 6 

127 ± 23 

3.43 ±0.18 




23.0 

24.0 

(1.10 ±0.16)- 

10"8 

87 ± 20 

3.20 ±0.16 




24.0 

25.0 

(6.91 ± 1.65) ■ 

10~ 7 

105 ± 28 

3.64 ± 0.44 




25.0 

26.0 

(9.28 ± 1.44) ■ 

10- 7 

139 ± 37 

3.35 ± 0.19 




26.0 

27.0 

(6.06 ± 1.34) ■ 

10- 7 

152 ± 44 

3.86 ± 0.48 




27.0 

28.0 

(7.34 ± 1.39) ■ 

10- 7 

126 ± 44 

3.27 ± 0.22 





a fixed 

b synchrotron energy flux 
c synchrotron vF v peak 
d e~ spectral index 


Table 2 

Results of the time-resolved spectral analysis of LAT data (100 MeV - 10 GeV) where TS is the value of the Test Statistic. 


T s tart [s] 

Tstop [s] 

TS 

7ph 

Te [erg s 1 cm 2 ] 

13.3 

23.7 

26.3 

— 2.9±0.7 

(7.4±3.5)xl0 “ 8 

23.7 

75.0 

26.8 

— 2.1±0.5 

( 1 . 6 ± 1 . 0 )xl 0 -8 

75.0 

237.1 

32.9 

—1.5±0.4 

(7.0±4.5)xl0 -9 

316.2 

724.8 

7.9 


< 2 . 0 x 10 -8 

3585.8 

6446.8 

1.8 


<7.7x 10 -10 

9307.8 

10000.0 

0.0 


<3.1xl0 -9 




Table 3 





The best-fit parameters of the E p 

evolution fit for each assumed value of no 

as well as the inferred r^. 


no [cm 3 ] 

To 

V 

9-3 

9 

r d [cm] 

1 

1125 . 9 );)®;® 

0 00" 1 " 0 ' 08 
u ‘ uu - 0.00 

0 17+ 0 - 01 
u ,±/ - 0.13 

i-25±g;l 

1.1 • 10 17 

10 

844 . 0 l )§; 3 

0 00 + °' 08 
u ' uu - 0.00 

0 1 7 + 0.02 
u,± ‘ - 0.14 

l - 26 ± S;ir 

6.0 ■ 10 1 ® 

100 

632 . 16 l 7 ; 7 

0 03“^ 0 05 
u - uo - 0.02 

0 1 7 + 0.05 

U - ±/ - 0.06 


3.4 ■ 10 1 ® 


for g < 1.5 in Equation [8] Since the obtained value of g 
via the fit to the E p evolution is not actually measuring 
the asymptotic behavior of the blast-wave evolution as 
is intended, we set g = 1.6 in Equation [8] corresponding 
to a nearly adiabatic blast-wave motivated by the low 
magnetic fields measured via t he s pectra and late time 
XRT observations (see Section [5~i] for full details). This 
choice of g does not affect the other parameters in the 
fit as g only modifies the late time evolution of E p (see 
Figure [l]). 


The analytic model given by Equation [7] and Equation 
[8] is consistent with the data up to a scale factor of 9 
which can easily be due to insufficient knowledge of all 
intrinsic parameters, on the simplified evolution law for 
the flux, and on the assumed value of g. It is pertinent 
to note that the parameters used in the P p (t ) prediction 
come from the fitted E p evolution alone. In addition, 
the predicted P p {t ) is independent of the observed flux 
values in the data. The close agreement between the 
predicted and measured P p {t) values is evidence that the 





















7 


7 io 53 

m 

be 

1 10 52 

-Q. 

CLh 

r-, 10 54 

I 

m 

S? 10 53 

<x> 


10 52 


Figure 6. The P p — E p plane of GRB 141028A (top) where P p 
is the i/F u flux calculated at E p . The predicted evolution ( green) 
is produced by substituting the recovered parameters of the E p 
evolution fit into Equation [T] The flux given by Equation [T] is 
scaled by a factor of 9 to match the data. At late times (pink 
region), the observed flux decays faster than the model predicts. 
The decay phase synchrotron HIC (blue) (bottom). The fitted value 
of the slope is £ = 1.33 =b 0.22. The decay phase is selected as the 
portion of the lightcurve that monotonically decreases with time 
i.e. from 10-27 s. 

observed flux at the peak and the observed vF v peak 
energy are linked in a manner predicted by the exter¬ 
nal shock model. We do note the deviation from the 
predicted curve at late times. This therefore provides 
strong evidence that these two independent quantities 
are linked in a manner predicted by the external shock 
model. 



E p [keV] 


4.3. The Luminosity - E p plane 

The commonly observed correlation in the decay phase 
of many GRB pulses between luminosity and E p in the 
form, 

(9) 


L oc E ^, 


sometimes r eferred to as the hardne ss-intensity correla¬ 
tion (HIC) ( |Golenets kii et al. 19831, is potentially an¬ 
other clue to the radiation mechanism responsible for 
the observed emission. It is quite simple to relate the 
luminosity and vF v peak energy of many common ra¬ 
diation mechanisms analytically. However, the spread 
in HIC power law slopes observed across many GRBs 
is difficult to explain with one mechanism (Borgonovo 
& Ryde 2001). However, single pulse GRBs lit with a 
synchrotron photon model have been found to have HIC 
p ower law slopes of ~ 1.5 without much spread (Burgess 


et al. 


2014b). 


W Uermer (2004), an analytic form for the HIC re- 
sulting from synchrotron emission is derived and pa¬ 
rameterized to account for various effects including ex¬ 
pansion geometry, magnetic flux-freezing, and radiative 
regime. Following this derivation, we can make the fol¬ 
lowing predictions for the evolution of the important 
quantities in determining the relationship between L and 
E p during the deceleration phase of a blast-wave. We 
have x oc tVCs+B, p oc i~g/{^g+1 ) ) £? p ^ TB"/F n , and 
L oc T 2 H 2 7 ^ in . Here 7 m i n is the minimum electron en¬ 
ergy in an assumed shock accelerated power law distri¬ 
bution (though the result is insensitive to the fact that 
we assume a power law) and B is the strength of the 
magnetic field. There are two regimes of cooling for elec¬ 


trons: fast and slow ( Sari et ah|l998 ). Each regime can 
be characterized by how 7 m i n evolves with time such that 


T 4 oc t~ 4 9/( 2 9+ 1 ) slow cooling 
(oT)” 1 oc i _2 9/( 2 9 +1 ) fast cooling 


Since L oc BTE p , we can write 


L oc 


p 2 slow cooling 
E p +9 fast cooling 


E, 


( 10 ) 


( 11 ) 


When fitting the L — E p data in the decay phase (10-27 
s) we find £ = 1.33 ± 0.22 (see Figure 6|. This is closer 
to the value expected for the slow-cooling regime which 
falls in line with the use of slow-cooling synchrotron to 
fit the time-resolved spectra. An a posteriori justification 
for the use of slow-coo ling synchrotron to fit the spectra 
is discussed in Section [5j] 


4.4. The photosphere 

We interpret the observed thermal component in the 
framework of photospheric emission, i.e., when the out¬ 
flow becomes transparent at a radius ~ 10 10 — 10 12 cm, 
i.e., much below r^. Because of highly relativistic mo¬ 
tion, the observed time delay between photons from the 
photosphere and those from the external shock is small. 
As the thermal component is not dominant, the errors 
in the calculation of outflow parameters derived in Pe’er 
et al. (2007) are large and we therefore did not use this 
formalism "herein. 

We note that the blackbody is not statistically signifi¬ 
cant (below 3-cr confidence level) in the first five seconds 
after trigger. However, photons from the photosphere are 
expected to arrive slightly before synchrotron photons. 
The delay between the photospheric photons and those 
at the peak of the external shock light curve can be esti¬ 
mated as (ra—r p h)/cFo ~ 9.5 s corresponding roughly td- 
Thus, we added five seconds to the duration of the photo¬ 
spheric emission. The identification of the photospheric 
emission can give important constraints. First, the du¬ 
ration of the photospheric component (Af p h) sets the 
width of the expanding outflow (W) at the photosphere 


(see Begu i et ak||2013 Begue & iyyani||2()l l| Vereshcha- 
giii||2U14p , 


Atph 
1 + z 


R ph W R 

2r 2 c + c 2r 2 c 


■ph 

fc 


( 12 ) 


where the first term accounts for the expansion time of 
the outflow up to the photosphere, the second term is 
the light-crossing time of the outflow and the last term 
is the angular timescale at the photosphere. For the ob¬ 
served duration at hand and for the Lorentz factor val¬ 
ues expected in GRB physics, the second term in the 
right-hand side of the equation dominates. It implies 
W ~ cAf ph /(l + z) ~ 1.35 • 10 11 cm. 

As the width W, the Lorentz factor T and the total 
energy of the blast wave are constrained, the state of the 
reverse shoc k can be a posterio r i stud ied. Following the 
discussion in Derrner & Menon (2009), we define 


54 


(JL) 4/a w 

V 300 / yv t 


E52 
A /2 V n 0 


V« 


( 13 ) 



























































such that a?NR = fy'd is the radius at which the reverse 
shock becomes relativistic. We find £ = 0.8 which im¬ 
plies that the reverse shock crosses the expanding outflow 
while it is only mildly relativistic. This justifies the use 
of Equation J3] in our treatment. We note the value of £ 
is the same lor each combination of T 0 and n 0 ( see Table 

0 >- 

5. DISCUSSION 

5.1. Compatibility with the Late-time X-ray/optical 
Emission 

Under the assumption of an external shock in the 
prompt phase, the afterglow is produced by the same 
process (that is the deceleration of a blast wave by the 
CBM), and we do not expect any break in the light curve, 
as already se en in several GRBs observed by Swift 
(Zhang et al. 2006). Unfortunately, the XRT observa¬ 
tions only cover a short time range beginning a few ~ 10 4 
s after the prompt emission (See Figure [7]). At most, we 
can say that XRT late time observations are compatible 
with the emission from a decelerating blast-wave in the 
semi-radiative regime, i.e., a region between 1.5 < g < 3. 

In addition, combining observations by XRT and 
Chandra, displayed in Figure [TJ a break can be iden¬ 
tified in the light curve. The exact time of the break is 
unconstrained but it ranges from 10 5 s < tb < 4.3 • 10 s s. 
Identifying the break as a jet break, the jet open ing angle 
can be estimated to be (Ghirlanda et al.||2004| 



Time [s] 


Table 4 

Values of eg from the afterglow and E p assuming different values 
of no as well as the location of the relativistic reverse shock (£). 


n 0 

To 

e B via v c 

e B via E p 

5 

1 

1125.9 

< 1.2 ■ 10“ 3 

~ 2.9 ■ 10" 9 

0.8 

10 

884.0 

< 2.7- 10“ 4 

~ 2.9 ■ 10“ 9 

0.8 

100 

632.2 

< 5.7- 10" 5 

~ 2.9 ■ 10 -9 

0.8 




Table 5 




Corrected Emitted Burst Energy (E co r) 



no [cm 3 ] 

tb (s) 

0j (degrees) 

E C or 

(erg) 

1 

1 • 10 s 

2.65 

1.0 ■ 

10 51 

4.3 ■ 10 5 

4.58 

2.9 ■ 

10 51 

10 

1 • 10 s 

3.53 

1.7- 

10 51 

4.3 ■ 10 5 

6.1 

5.2 ■ 

10 51 

100 

1 • 10 s 

4.71 

3.1 ■ 

10 51 

4.3 ■ 10 5 

8.15 

9.2 ■ 

10 51 


50 ks after the trigger. The spectrum is consistent with 
a single power-law of index j3 = 1.29±0.07. Therefore, we 
can deduce that the cooling frequency is above the XRT 
frequency: v c > i/ X r T - From the expansion of a blast- 
wave in a constant d e nsity C BM, u c can be estimated as 
(Panaitescu & Kumarjj2000a) 

z'c = 3.7 • 10 16 E~ 3 1/2 n^(Y + 1 )- 2 e-%T~ 1/2 Hz, (15) 

where cb parameterizes the magnetic field in the shocked 
ISM, Xd is the time in days and Y is the Compton pa¬ 
rameter, that we chose to be zero for simplicity. This 
leads to an upper limit ee < 2.3 • 10 -5 taking T d = 1 day 
and assuming no = 100 and U 53 = 100 which implies 
T 0 = 632 (see however, Table E] for values corresponding 
to different parameter choices JT 

Also, this value can be cross-checked by considering 
the peak energy in the first seconds of the prompt phase, 
which can be evaluated as: 


Figure 7. The X-ray and optical afterglow of GRB 1^1028A. 
A power law fit best describes the optical data while a broken 
power law fit, indicative of a jet break, best describes the X-ray 
data. The yellow shaded region indicates the 90% uncertainty in 
the break time due to the sparse sampling of data prior to the 
Chandra follow-up observation. 


dj = 0.161 


tb 

1 +z 


3/S 


77 7 n 0 


E, 


7 ,iso,52 


i/s 


(14) 


where tb is the break time in days and ry 7 is the radiative 
efficiency taken to be 0.1 corresponding to our assump¬ 
tion of Eq = 10 55 erg. Thus, the opening angle ranges 
from 2.65° < 9j < 8.15°, whe n con sidering the fitted pa¬ 
rameters obtained in Section |4.1| Therefore, the radia¬ 
tive energy budget of the burst is considerably reduced 
to few ~ 10 51 erg (see Table [ 5 J. 

Finally, additional constraints can be obtained from 
simultaneous optical and X-ray observations at around 


E p ^ 


J32TTeBm p c 2 n 0 T 2 2 
-R- r7mi 

-Drr it. 


(16) 


where B cr it = 4.414 • 10 13 G is the critical magnetic field, 
k parameterizes the minimum Lorentz factor 7 m i n of the 
accelerated electrons such that 7 m j n = «T( m p/m e ), m e is 
the mass of an electron, and E p is given in units of the 
mass of an electron. Using E p ~ 1.3 MeV and no = 100 
as obtained from the data, and assuming n = 1 gives 
cb ~ 3.0-10 -9 (see also Table I| . This value is consistent 
with the upper limit obtained from the late afterglow^ 
Finally, the cooling time of an electron with Lorentz 
factor y m i n can be compared to the expansion time of 
the blast-wave at the observed luminosity peak which 


9 However we note that this value is very sensitive to k which 
is fairly unconstrained: as an example, with k = 0.1, it becomes 
£b ~ 2.8 ■ 10“ 5 , incompatible with the upper limit obtained from 
the afterglow. 
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is on the order of the id- We find that they are com¬ 
parable, leading to efficient energy extraction from the 
electrons, without drastically changing the electron dis¬ 
tribution function, *.e., creating an additional power law 
over several orders of magnitude at low-energy charac¬ 
teristic of fast-cooled (or completely cooled) electrons. 


5.2. LAT late time GeV emission 

Fermi LAT observed a GeV component over the du¬ 
ration of 10 — 10 3 s. Unfortunately, there were no simul¬ 
taneous observations at other wavelengths from which a 
broadband spectrum could be obtained. Nevertheless, 
we can use the analytic flux evolution derived in Der- 


mer et ah (1999) to check if the observed LAT fluxes are 
compatible with the model. We must first address a few 
caveats: the analytic model assumes no synchrotron self- 
Compton (SSC) emission and the fluxes are sensitive to 
our lack of knowledge about the total burst energy (which 
we assume to be 10 55 erg for the purpose of calculation) 
as well as the degeneracies in the recovered parameters 
from the E p evolution fit. 

In the past years, a debate concerning the origin of 
the late time (few seconds after the trigger) high-energy 
LAT emission has taken pl ace, o ppos ing a synch rotron 
mechanism (|Kumar fe Barnio l Duran 20091 2010), SSC 
dPanai tescu &E Kumar 2000 b|, and pair-loading of the 
CBM (Beloborodov 2005) as the possible candidates for 
this emission. First, we can assess the SSC component, 
and show that with the parameters at hand from the 
prompt emission the SSC flux in the LAT bandpass is 
much less than the synchrotron flux. We can constrain 
the SSC flux in the LAT band p ass by following the dis- 
Beniamini & Piranf ([2 010 (see Appendix |B| 


cussion m , _ 

for further details). Via Equation B12| the ratio of S 


flux to synchrotron flux in the LAT bandpass (F ssc ) can 
be computed by substituting values derived in Section 
5.1 to yield F ssc ~ 10 -8 . The numerical result shows 
that the SSC flux in the LAT band is much less than 
the synchrotron flux. However, the computation relies 
on the assumption that electrons are fast cooled either 
by synchrotron or by inverse Compton of the produced 
synchrotron photons, that is to say that all the energy is 
emitted by one or the other mechanism. It comes with 
two consequences: 


1. th e Com pton para met er is an upper limit as shown 


Sari et al. (19961. Therefore the result of 


Eq.(B12) is only an upper limit, 


2 . as a result of the assumption of fast cooling, 
the electron distribution function extends to small 
Lorentz factors. These electrons are not in the 
Klein-Nishina regime and can efficiently upscat- 
tered synchrotron photon in the LAT band. How¬ 
ever, we found that electrons with Lorentz factor 
7 min are not strongly cooled over a dynamical time. 
Therefore the electron distribution function does 
not extend substantially below 'ymin, reducing the 
inverse Compton flux. 


We now consider that the GeV emission is the result 
of synchrotron emission from the same forward-shock re¬ 
sponsible for the lower energy prompt phase. Using the 
analytic estimate for the blast-wave evolution, we esti¬ 


mate the expected flux of the synchrotron emitting blast- 
wave and compare it to the integrated (100 MeV - 10 
GeV) L AT ener gy flux (Fe) by integrating Equation 1 
from Dernier et al. (1999) across the LAT bandpass (see 
Figure [0. Because of our lack of knowledge about the 
intrinsic total energy of the blast-wave and the degenera¬ 
cies in the fit parameters (g, n o), we vary the parameters 
across a broad range and find that the model gives con¬ 
sistent limits within an order of magnitude. We also 
note that this model is simple and will not be as accu¬ 
rate as a full numerical solution. Additionally, using the 
power-law electron index ( p ) found from GBM+LLE syn¬ 
chrotron fits, we compute the synchrotron photon indices 
via the transformation y p h = — (p + l)/2 and compare 
them with the photon indices found in the LAT alone 
(Figure [0. There is a clear evolution of the synchrotron 
photon index increasing that smoothly transitions to the 
indices observed by the LAT. 



Figure 8. The LAT GeV integrated energy fluxes (100 MeV - 10 
GeV) are com pared with a range o f expected fluxes from the ana¬ 
lytic model of |Dermer et al. (1999]) that includes only synchrotron 
emission. Suitable agreement can be found for a valid range of 
assumed blast-wave parameters. 



Time [s] 


Figure 9. The photon indices from the GBM+LLE fits and LAT 
fits. A clear evolution is seen into the LAT band pass with the 
indices all in line with what could be expected from synchrotron- 
emitting high-energy electrons accelerated into a power law 
momentum-distribution. 


Therefore, we conclude that the high energy emission 
from the LAT is compatible with synchrotron emission 
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alone, and does not r equi re an a dditional me chanism. It 


was proposed by Kumar & Barniol Duran (2009 
in which the prompt IVleV phase is assumed to 


2010 ), 
De the 


result of an unspecified mechanism, that the GeV LAT 
emission arises from a separate exernal shock with a syn¬ 
chrotron peak on the order of a few 100 MeV. Conditions 
on the density and the magnetic field are then derived 
such that this external shock is not more luminous than 
the prompt emission in the sub-MeV band. However, for 
GRB 141028A we show that the entire emission can be 
explained by synchrotron emission with an MeV peak en¬ 
ergy from a single forward shock without fine-tuning of 
the parameters, i.e., both the MeV and the delayed GeV 
can be explained self-consistently by an external shock. 


5.3. Summary of Prompt Phase 

In this section, we concisely summarize the novel ap¬ 
proach we have employed to test the external shock 
model in the prompt phase. Almost as important as 
the observed E p evolution is the fact that we can fit 
the spectrum with a physical slow-cooling synchrotron 
model. Even though it is possible to fit a physical syn¬ 
chrotron model to the spectra (as opposed to an empir¬ 
ical Band function from which the comparison to phys¬ 
ical models must be inferred and c an lead to problems 
as shown in Burgess et al. (2014a)), it must be shown 
that the photon model is sound, i.e, that the dynamical 
evolution can lead to slow-cooling synchrotron emission. 
In Section [O] we confirmed that both the afterglow and 
prompt phase are consistent with slow-cooling emission. 

The evoluti on o f E u in ti me can be fit with the analytic 
prediction of Dermer et al. (1999). While the blast-wave 
appears to be transitioning to the asymptotic limits, (g < 
1.5) we tentatively conclude that it is evolving towards 
an adiabatic regime. Full temporal coverage of the flux 
and spectra through the late time evolution will help to 
resolve this issue. Unfortunately, few single pulsed GRBs 
have been observed with such a temporal and wavelength 
coverage. This unique opportunity presents itself only 
due to the current multi-mission capabilities and begs 
for continued multi-wavelength coverage. 

If we take the fitted parameters of the E p evolution fit, 
we can predict the evolution of the P p — E p plane which is 
consistent with the data in both the rise and decay phase 
of the pulse. The HIC of the prompt decay phase of the 
synchrotron luminosity is consistent with what would be 
expected from slow-cooling synchrotron in a decelerating 
external shock. 

We use the observation of the photospheric component 
to estimate both the width of the blast-wave as well as 
when the relativistic reverse shock crosses the forward 
shock. From our calculations, the reverse shock will have 
very little impact on the evolution of the forward shock 
and we therefore neglect it. 


6 . CONCLUSION 

In this article, we have shown that the prompt phase 
of GRB 141028A is consistent with originating from an 
external shock that emits synchrotron radiation from an 
electron distribution that did not have enough time to 
cool completely (slow-cooled). Combining the results of 
the E p evolution fits with the other clues from the data, 
an external shock origin of GRB 141028A is a very likely 
scenario. Not only do we find that the prompt emission is 


explained by synchrotron emission from a forward shock, 
but we find clues for a late time high-energy emission 
resulting from the extension of this synchrotron emis¬ 
sion rather than an SSC component. We want to stress 
that in this scenario the delineation between “prompt” 
and “afterglow” is superfluous as the early and late time 
emission both originate from the same mechanism. We 
are currently applying this analysis to other long, single 
pulse GRBs and find that they too have their vF v peak 
and P p {t) which evolved consistently with the predictions 
of an external shock. These results will be presented in 
a forthcoming publication. 

We propose that external shocks are still a viable can¬ 
didate to explain the prompt dynamics of GRBs. How¬ 
ever, the most likely scenario is that they are a subset of 
multiple dynamical frameworks including internal shocks 
that further subdivide into categories based on opacity 
and/or magnetic content of the outfl ow. We no te t hat 
this is not an entirely new idea. : Panaitescu & Meszaros] 
(19981 numerically investigated the parameter space of 
the external shock model and proposed that smooth, 
single-pulsed (and possibly multi-pulsed) GRBs can be¬ 
long to a subclass of GRBs that are the results of exter¬ 
nal shocks. We have now shown this quantitatively by 
applying the model to the data. 

The analysis of individual GRBs is crucial to identi¬ 
fying what is likely a host of different emission mech¬ 
anisms. For example, GRB 090902B exhibits dynamics 
and spectra that are consis tent with subphotospheric dis¬ 
sipation (R yde et al.[ i2010). The fact that both external 
shocks and photospheric emission are observationally vi¬ 
able candidates for explaining GRB emission implies that 
we should use caution when trying to apply one mecha¬ 
nism as an explanation for all GRBs based of properties 
from catalogs and should instead focus on how to relate 
the different emission mechanisms into a unified frame¬ 
work. 
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APPENDIX 


LAT ANALYSIS 

The model used in the likelihood fit is composed of the Galactic diffuse emission produced by cosmic-ray interaction 
with gas and radiation fields and the isotropic diffuse emission . In additi on, w e add all the point sources in the ROI with 
spectral models and parameters from the 2FGL catalog (Nola n~et al.|||2012 ). While the normalization of the Galactic 
template and the spectral parameters of all the 2FGL sources are frozen to their nominal values, the normalization 
of the isotropic template is left free to vary in order to absorb statistical fluctuations. The GRB location is fixed, 
and its spectrum is described by a power law d N/dE oc -E 7ph with y p h the photon index (note that typically 7 p h <0). 
Following the time-resolved analysis described in Ackermann et al. (2013) we split the LAT data in 48 log-spaced time 
bins from 0.01 to 10000 seconds after the trigger. The logarithmically-spaced binning provides constant-fluence bins 
when applied to a signal that decreases approximately as 1/time, such as the extended GRB emissions observed by the 
LAT. We first merge consecutive time bins in order to have at least 5 counts per bin (corresponding to the number of 
free parameters in the likelihood model plus 2), and we then perform likelihood analysis. We estimate the significance 
of the GRB source by evaluating the “Test Statistic” (TS) equal to twice the logarithm of the ratio of the maximum 
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likelihood value produced with a model including the GRB over the maximum likelihood value of the null hypothesis, 
i.e., a model that does not include the GRB. The probability distribution function (PDF) of the TS under the null 
hypothesis is given by the probability that a measured signal is compatible with statistical fluctuations. The PDF in 


from Wilks’ 

theorem i 

Wilks 

19381 Protassov et al. 

2002 

). However, it has been verified by dedicated Monte Carlo 

simulations ( 

Mattox e 

EaT . 1 

961) that the cumulative PDF ol the TS in the null hypothesis (i.e., integral ol the TS 
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degrees of freedom associated with the GRB. The factor of 1/2 in front of the TS PDF formula results from allowing 
only positive source fluxes. 

If the resulting TS value is lower than an arbitrary threshold (TS < 10) we merge the corresponding time bin with 
the next one, and we repeat the likelihood analysis. This step is iterated until one of two conditions is satisfied: 1) we 
reach the end of a GTI before reaching TS = 10, in which case we compute the value of the 95% CL upper limit (UL) 
for the flux evaluated using a photon index of —2; 2) we reach TS > 10, in which case we evaluate the best-fit values 
of the flux and the spectral index along with their ler errors. 


HIGH-ENERGY SSC FRACTION 


Following Beniamini & Piran (2013), the SSC and Klein-Nishina peak frequencies are defined respectively as 


Cnin — E p h 
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where h is Planck’s constant. The total SSC flux is linked to the total synchrotron flux at the peak (F/ 
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where Y is the Compton parameter and, 
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In the derivation of the Y via Sari et al. (1996), the Klein-Nishina correction of the Compton cross-section (<7 KN oc 
ln(2a:)/a; where x is the photon energy in the rest frame of the electron in m e c 2 units) was ignored. Here, in order 
to estimate Y we include this correction by considering electrons with Lorentz factor 7 m i n and their corresponding 
synchrotron photons. Therefore, we write 
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The total upscattered flux in the LAT bandpass can then be written as 
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Here, a refers to the low-energy photon index of the Band function which is -2/3 in our case and y nlM is the maximum 
synchrotron frequency corresponding to the maximum Lorentz factor of the accelerated electrons (de Jager et al. 1996): 
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In order to estimate the synchrotron flux in the LAT band pass, we write 
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where s is the high-energy electron power law index which we take s = 2.5 for simplicity. Then we can write the 
fraction of SSC to synchrotron flux in the LAT bandpass as 


qsc — 


PLAT 
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T SSC 


f dv ( 

J LAT y ^min J 




This can be further simplified via Equations B4 and |B7| to 
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This is in fact an upper limit on the fraction of SSC due to the assumption that electrons radiate all their energy in 
a dynamical time either by synchrotron and/or inverse-Compton and that there is a significant fraction of electrons 
below 7 m i n . If either assumption fails (as is true in our analysis) the fraction of SSC emission will be suppressed. 






